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Abstract 

A recent proposal that gravity theory is an emergent phenomenon also entails the 
possibility of photon decay near the Schwarzschild event horizon. We present a possible 
mechanism for such decay, which utilizes a dimensional reduction near the horizon. 
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In a recent paper the question is raised whether classical General Relativity can 
exist as an emergent phenomenon — as the low-energy limit of an underlying quantum 
system. In this view, the singularity at the Schwarzschild event horizon represents a 
failure of the effective description owing to the divergence of a characteristic coherence 
length. We shall introduce an additional element to this model: a natural change in the 
dimensionality of virtual particle loop integrals of the quantum system near the horizon. 
This change will have important implications. In particular, it may cause photons to 
decay when they near the event horizon of a black hole, as suggested in 

In the Schwarzschild metric, the line element is 



r 



(The velocity of light is set to one.) Classically, this metric concentrates the motion in 
the radial direction near r = r^. This can be seen by examining the spatial line element 



dt'={i-!^ 

r 
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dr^ + !-!!.] r^dQ^ 



(2) 



For r near rs, the angular variables are suppressed, and motion is confined to the two- 
dimensional t-r subspace. 

More specifically, consider the geodesic equation, + F^^x'^i;^ = 0, where x'^ is the 
derivative of the position with respect to an affine parameter. 

The temporal and angular equations, 

i+—^-^—ir = 0, (3) 
r[r — rs) 

9 + 2'^9-^sm29<p'^ = 0, (4) 

4)+-r(p + 2cotee<j) = 0, (5) 
r 

can be readily integrated to give 



■ r-rs , . 

t = r, (6) 



r 

sin^ 9 (f) = m, (7) 

ir'er + ^, = I'- (8) 
sm 9 

r, m, and P are integration constants. Near the Schwarzschild horizon we may allow all 
particles to be massless and take the geodesies to be null. Then the constants will diverge, 
but we retain their ratios, which remain finite. From the 9 equation, it is clear that for 
9 = ^ and 9 = 0,9 = also, so the motion remains in the equatorial plane, with = P. 
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Using the other three equations, the radial equation 

+ (i - + - t) " ^''^ + - + " = (8) 

becomes 

fH ^ ^ ^ ^ = 0. (10) 

2r(r — Ts) 2r{r — rg) 

In this form, the equation may be integrated to give 



±rWl- (1--^^, fill 



^ . 2 ' 



where X = 1/t and a possible integration constant is set to zero for null geodesies. 

Dividing the (p and r equations, (|^) and (PTD, by the t equation (^, gives at 6' = |, 
m = /, 



r5\ A 

(it V r J r 



^--)-2^ (12) 



Consider a photon trajectory just outside the event horizon, at radius r = r5 + A. For 
A -C r5, the time-evolution equations (|12D and ([T3|) become 



A 



37 = 13^- (15) 



dt 

The equations, with the lower sign in (p!^, may be integrated, yielding 

cosh 



A = (16) 



2rs 

= ^tanh— . (17) 

Equations (16) and ( p!?]) show that if ^ is sufficiently small (^ tt), then the radial motion 
will be much more rapid than the angular motion. This is the regime we are interested 
in, where the concentration of motion in the radial direction reduces some aspects of the 
system to 1+1 dimensions. 

We shall consider a much more profound reduction in the dimensionality of the system. 
Taking local coordinates r, x, and y, the spatial line element is 

df = -grrdP + dx^ + dy"^. (18) 
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In the model the singularity at r = represents a real physical effect, not merely a 
coordinate artifact, so these are very natural linear coordinates. If we suppose that p^, 
p^, and are cut off at the same scale in a loop integral, = — p^p* will be dominated 
by —grr{p^Y- (This supposition strongly breaks general covariance, of course.) This is 
the kind of situation we would like to analyze. However, a sharp momentum cutoff is 
not gauge invariant; to study the photon self-energy, we shall translate this idea into 
gauge-invariant language, using dimensional regularization. 

In the renormalization prescription described above, one of the spatial dimensions 
provides the dominant contribution to p^. In the language of dimensional regularization, 
this can be seen as a reduction in the effective dimensionality d of the momentum integral 
to (i < 4. 

We shall find the effective dimensionality by examining the volume element of this 
system, because the momentum cutoff in a given direction and the volume contribution 
of that direction are closely related. To see this, consider for the moment a theory that 
is regulated by a lattice at short distances. The volume of a region counts the number of 
lattice points in that region. The lattice spacing Oj in a given direction governs the density 
of lattice points along that axis, so a length L in the Xj-direction contributes an amount 
to the volume, where a is a scaling constant independent of direction. The lattice 
spacing also corresponds directly to the momentum cutoff in that direction, p"^"-^ = 2L, 
So the dependence of the volume on a given dimension and the momentum cutoff in that 
direction are intimately linked. 

So we look for an expression for the effective dimensionality (to be used in dimensional 
regularization) in terms of the volume element. In Schwarzschild space-time, the spatial 
volume element is 

dV = ,/g dr de d(j). (19) 

Here, g = —grrf'^ sin^ 6 is the determinant of the spatial metric —gij- The classical geodesic 
problem suggests that the radial direction should always contribute one effective dimen- 
sion, while the angular directions may contribute less than one. We wish to determine the 
effective dimensionality by an integral over an effective volume. To a achieve a reasonable 
result, the effective volume is defined by rescaling a factor of \/—g^ from each direction. 
So the new volume element reads 

T"^ sin 

dV = dr de d(t) = gttdVE, (20) 

9rr 

where dVE is a Euclidean volume element. 

From dV, we need to find an expression for the number of effective dimensions. This 
expression should have several properties. The dimension corresponding to a volume 
element dxi dx2 ■ ■ ■ dxn should be n. (with dVE corresponding to 3 dimensions). So the 
dimension function should be additive where the volume element is multiplicative; this 
is the fundamental property of a logarithm. So a natural choice for the effective spatial 
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dimension dg is 

InA • 

Evaluating this gives us 

rA 



(21) 



= Wo^ + 2^ = 3 + 2^^. (22) 
InA InA InA ^ ^ 

Equation ( p2l) has some problems. The most striking one is that A appears to be a 



dimensional quantity, which would make ambiguous. The obvious solution is that 

the coordinates must be nondimensionalized, to make A dimensionless. However, it is not 
at all obvious how to nondimensionalize the coordinates. Fortunately, we do not need 
to deal with that question directly. Regardless of the coordinates' dimensions, dVs will 
always contribute three dimensions to dg- To analyze the last term in (p^), we introduce 
the natural condition that the angular factors in the volume element can not contribute 
any less than zero dimensions each. That is. 



In ,/g77 , , 

> -1. 23 

InA 

Since Qu < 1 and the logarithm is strictly increasing, this condition may be rewritten as 

< VaTf (24) 

In classical General Relativity, gtt{rs) = 0, but in the model [jl|, gu drops to a nonzero 
minimum value 

3tt ~ 1 ~ — > (25) 

rs + rs 

where 5 is a small length that characterizes how the classical Schwarzschild singularity 
is cut off by the underlying quantum system. We expect 6 to be related to the Planck 



length. Thus we estimate A ~ yrs/S- Close to the event horizon, at r = + A, where 
5 ^ A ^ r5 and gu ~ ^/fs, the effective dimension of the system is 

So we will look at the problem of Quantum Electrodynamics in (1 + e) + 1 dimensions, 
where 

We shall only consider the contribution to the photon self-energy from massless par- 
ticles. This should be a good approximation near the event horizon. Photons coming 
in from spatial infinity are highly blueshifted at r = + A, so the momentum scale in 
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the photon propagator is large compared to any invariant momentum scale (such as the 
electron mass). The same argument may also be phrased in different terms. The energy 
of a comoving electron of mass rrie is y^rrie, so near the event horizon, the apparent 
electron mass becomes small. So it is reasonable to consider massless particles. 

The one-loop photon self-energy due to a single species of charged massless fermions 

is 

tU'^-'iq) = 2trI(gV^ - g^) f dx [ (28) 

Jo Jk [k'^ + q''x[l — x)\^ 

Here, I is the unit matrix in spinor space, and the the Minkowski metric is denoted by 
rj^'^, to avoid confusion with the GR metric g^". Since d ^ 4 as r —>■ oo, the Dirac 
matrices should be four- dimensional. Equation (ESp, derived in the appendix, gives, upon 
a (i- dimensional /c-integration. 

There is a subtlety in the use of (^). In dimensional regularization, it is usual to 
reduce all aspects of the problem to d dimensions. In our case, only the loop integral 
is c?- dimensional. There are still four Dirac matrices, and the photon remains a four- 
component field. However, so long as the external photon momentum q lies in the d- 
dimensional subspace, equation (BO) remains valid, and the metric rj'^'' is (i-dimensional. 



The d = A and d = 2 cases of equation (^9]) are well understood. Since we have 



c? = 2 + e, we expand around d = 2. Evaluating equation ( poD with this value of d, we get 

2p2 1 f7ri/2-e/2 /o 1 ^ 



The bracketed term in (pOD is unity at e = 0. It is purely real, so will only contribute 
higher-order corrections to the real and imaginary parts of the self-energy. 

At d = 2, we get Schwinger's well-known result that the photon becomes massive 0. 
The self-energy is 

^n'^^(^)l.=2 = ^v(^''~l?)- ^^^^ 

(This differs by a factor of two from the usual result, because here we have used four- 
dimensional Dirac matrices.) The residue of the pole at = gives the photon mass 



For d = 2 + e, the result is only slightly different. Instead of having we have 
(^_q2y-^/2 i which we expand about e = to get 

1 ] 



.q2y-e/2 _g2 



l + -ln(V) 



(32) 
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As in equation (|3TD, there is a pole at = 0. However, the residue is shifted by the 
second term in (p2[). To zeroth order in e, the pole in the propagator is shifted to To 
first order, the pole is further shifted to the value of [1 + | ln(— g^)] evaluated at the new 
pole location. This shifts the pole to 




(33) 



Equation (|33|) is correct to lowest order in e in both its real and imaginary parts. The 
sign of the imaginary part has been chosen so that photons decay rather than appear. 

We must now turn to the question of how to interpret equation (^31). By expanding 
around = 2, we have introduced a number of two-dimensional conventions. The 
appearing in (^) is the two-dimensional value of the electromagnetic coupling. In two 
dimensions, e has mass dimension one, so — is indeed a mass squared. We must relate 
the e in ( ^3]) (which we will henceforth refer to as 62) to the four- dimensional electron 
charge 64. 

We may relate the two charges by comparing the actions in two and four dimensions. 
In four dimensions, the electromagnetic Lagrangian density is 



£4 



Then the action is 



5*4 = y dtr^dr d^l Ci. 



(34) 
(35) 



The action S2 derived from the two-dimensional Lagrangian C2 should the same as S4, 
up to a constant factor C. So we have 



/ 



dtdrC2 = —C 



dtPdr dn^F^^F^"". 



(36) 



We must perform the angular integrals on the right-hand side of (^) to determine £2- 
This means doing an integral over the submanifold orthogonal to the t-r subspace. This 
orthogonal submanifold is a sphere, and the integral over it will depend upon the radius 
at which the integral is evaluated. We are interested in in radii r rs (which is the only 
region where the integration over angles is justified). If F^^ is spherically symmetric, the 
angular integral gives 47r, and we can set r = to get 



j dtdrC2 = -4nrlC J dtdr^F^^F>'\ 



(37) 
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We can now read off the value of £2, 



= (38) 



so that the two-dimensional charge is 

1 



- Anrlc'"'- ^^^^ 



The constant C includes differences in how the field operators are normalized in two and 
four dimensions. So although the precise numerical relation between 62 and 64 has not 
been determined, the dependence of 62 on is unambiguous. 

In deriving equation (|5UD, we assumed that the field configuration was spherically 
symmetric. We can also evaluate the angular integral for more general field configurations, 
although this adds additional ambiguities. If the field F'^" is in an / > 0, m = multipole 
state, the integral becomes 

J dtdrC2 = -C J dtr'^drdnPi{cos9f^^Ff,^F'"'. (40) 

(The additional angular fields coming from derivatives of Pi{cos6) are suppressed by 
^ and have been dropped.) Since the maximum value of Pi{cos6) is Pi{l) = 1, the F^'^ 
appearing in equation (|40|) is the maximum value of the field over all angles. It is consistent 
with our earlier identification of the two- and four- dimensional fields to identify this F^" 
with the field appearing in £2 although other identifications could also be consistent). 
Evaluating the integral then gives us 

A similar calculation may be done for m 7^ 0, but the result (with these conventions) 
depends on m explicitly, not merely on /. Despite this problem, (^l]) remains a good 
candidate for an m-independent multipole field mass. 

We must also address the question of how to interpret the imaginary part of (0). 
To help with the interpretation, we shall use an analogy to a much simpler dimensional 
reduction problem — the electromagnetic field in a rectangular waveguide This simple 
problem in classical electrodynamics has many similarities to the QED problem under 
consideration. 

Consider a rectangular waveguide with metal walls. The waveguide has dimensions a 
in the x-direction and b in the y-direction. (We will presume that a and b are comparable 
in magnitude.) The waves propagate freely in the 2;-direction. The boundary conditions 
on this system restrict the wavevector of the electromagnetic field in the interior to be 

k=^x+— ^y + M. (42) 
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The numbers rix and Uy are positive integers; at least one of Ux and Uy must be nonvan- 
ishing for fields to exist. The frequency a; = |k| satisfies 



Since propagation only occurs along the z-axis, it is natural to look at this system in the t- 
z subspace, where the wavevector is simply kz- Then (^) looks like the energy-momentum 

relation for a relativistic particle of mass ml^g = vr^(^ + -^). 

So in 1 + 1 dimensions, a photon in a waveguide acquires an effective mass. The 
scale of this mass is a~^, where a is the characteristic size of the system in the neglected 
dimensions. This is the same scaling we found previously. In equation (|33D , the scale of 
the photon mass was r^^, and rs is the length scale of the event horizon in the angular 
directions. According to the black hole system actually has a whole hierarchy of 
photon masses. The waveguide also exhibits this property; different and Uy values give 
different values of m^^. (These results are similar to those found in Kaluza-Klein theories, 
although the higher modes are not strongly suppressed here.) 

The waveguide system also exhibits another important property — decay. Through 
interactions in the x- and ?/- dimensions, a photon can disappear from the interior of the 
waveguide. This can occur in a variety of ways, depending on the regime. We mention 
the two simplest regimes and discuss the interpretation of decays in these regimes. At low 
frequencies, u <^ u, where u is the collision frequency for electrons in the metal walls, the 
magnetic field drives surface currents which dissipate energy through resistive heating. 
This leads to a low-frequency energy loss 

where a is the conductivity of the metal walls. The behavior is different at higher frequen- 
cies, u <^ Lu < Up, where Up is the plasma frequency, related to the electron density Ue by 
Up = In this regime, the electromagnetic field is exponentially damped in the walls, 
but photons can tunnel through the walls and escape from the waveguide. However, the 
decay rate does not have a simple dependence on a and Up. 

The decay rate in the waveguide depends strongly on the regime, but in each regime, 
the decay rate depends primarily on the length a and some other length parameter. In 
the regimes outlined above, the length parameters are provided by a and Up. In the black 
hole model, the imaginary part of is governed by e, which depends on the inverse mass 
scale rs, as well as on the lengths 6 and A. 

The two decay regimes outlined above possess very different decay processes. The 
first regime is dissipative. Through interactions, photons in the waveguide decay into 
something else — thermal excitations of the metal walls. In the second regime, the photon 
does not actually decay. Instead, it escapes from the waveguide and the corresponding 1 + 
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1- dimensional subspace. Despite their differences, these processes would each contribute 
an imaginary part to wi^^. 

Analogously, the imaginary term in equation (|33|) could have two different origins. The 
photons could be decaying into other particles, as is suggested in Alternatively, the 
effective decay rate in (|33D could correspond to photons being scattered out of the 1 + 1 
dimensional t-r subspace into states with large angular momenta. Our simple treatment 
does not allow us to distinguish between the two. However, either process would be 
novel — an effect caused by the finite minimum of gu and controlled in magnitude by 6. 

Although it is not directly relevant to the problem at hand, it is worth mentioning 
one further suggestive aspect of the waveguide analogy. In 2 + 1 dimensions, the photon 
self-energy ( ^Qf ) is still singular at = 0, but that singularity is weaker than a simple 
pole in |Q. This introduces ambiguity as to whether or not the photon has mass in 
2 + 1 dimensions. The waveguide analogy of 2 + 1-dimensional QED is the propagation 
of the electromagnetic field between two parallel plates of separation a in the x-direction. 
However, in this case, the wavevector in the x-direction, is allowed to vanish. So a 
photon in this system may or may not behave as if it were massive. 
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Appendix: Photon self-energy integral 

To derive ( ^8]) in a dimension-independent manner, we begin with the integral 

Introducing a Feynman parameter x, setting = — x(l — x)q'^, shifting the integration 
variable k ^ k — xq, and dropping all odd-k terms leaves 

Evaluating the traces for D-dimensional Dirac matrices and a c?- dimensional integration 
over k gives 

= tr (^-7 + ^7^^^') 
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- - 1 ) Dt]'"'^ 



tr(7^ #7^ if) 



(47) 
(48) 



The integrand with the k"^ from (|47|) may be transformed to resemble the rest of the 
integrand according to 



L 



k' 



1 



1 



1 + 



M2 



k {k^-M^y \d 



(49) 



The first term of (|49|) may be Wick-rotated to Euclidean space and evaluated using 
integration by parts, to give 



fe P -M2 \d 



2 

- 1 



dQnk'^-^dk (2-d^ 



-A;2 _ M2 
dVtddk A;2 



A;2 + M2 



{d-2)k 



d-3 



n 



dVlddk k 

M2 



d-2 



d 



dk k^ + M2 



Adding this to the other term in ( ^9]) gives 

P /2 _ ^ 



fc (P - M2)2 \d 



d Jk (F + M^y 



M2 



(50) 



(51) 



Combining (ETf) with the integrals having numerator 
equation (|28|). 



k {k^-M^y 

gives the self-energy result 
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